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ABSTRACT 

We present Hubble Space Telescope observations of active asteroid 
313P/Gibbs (formerly P/2014 S4) taken over the five month interval from 2014 
October to 2015 March. This object has been recurrently active near perihe¬ 
lion (at 2.4 AU) in two different orbits, a property that is naturally explained 
by the sublimation of near surface ice but which is difficult to reconcile with 
other activity mechanisms. We find that the mass loss peaks near 1 kg s“^ 
in October and then declines over the subsequent months by about a factor of 
five, at nearly constant heliocentric distance. This decrease is too large to be 
caused by the change in heliocentric distance during the period of observation. 
However, it is consistent with sublimation from an ice patch shadowed by local 
topography, for example in a pit like those observed on the nuclei of short-period 
comet 67P/Churyumov-Gerasimenko. While no unique interpretation is possible, 
a simple self shadowing model shows that sublimation from a pit with depth to 
diameter ratio near 1/2 matches the observed rate of decline of the activity, while 
deeper and shallower pits do not. We estimate the nucleus radius to be 700T100 
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m (geometric albedo 0.05 assumed). Measurements of the spatial distribution of 
the dust were obtained from different viewing geometries. They show that dust 
was ejected continuously not impulsively, that the effective particle size is large, 

~50 fim, and that the ejection speed is ~2.5 m s“h The total dust mass ejected 
is ~10^ kg, corresponding to ~10“^ of the nucleus mass. The observations are 
consistent with partially shadowed sublimation from ~10^ m^ of ice, correspond¬ 
ing to ~0.2% of the nucleus surface. For ice to survive in 313P for billion-year 
timescales requires that the duty cycle for sublimation be <10“^. 

Subject headings: minor planets, asteroids: general — minor planets, asteroids: 
individual (313P/Gibbs (2014 S4)) — comets: general 

1. INTRODUCTION 

313P/Gibbs (formerly P/2014 S4), hereafter called “313P”, was discovered on UT 2014 
September 24 as a product of the on-going Gatalina Sky Survey (Gibbs 2014). Although 
313P is cometary in appearance, its orbit lies within the asteroid belt, having a semimajor 
axis 3.156 AU, eccentricity 0.242 and inclination 11.0°. Perihelion and aphelion occur at 
2.392 AU and 3.920 AU, respectively. The Tisserand parameter with respect to Jupiter, Tj 
= 3.13, is distinct from those of Kuiper belt and Oort cloud comets, which have 2 < Tj < 
3 and Tj < 2, respectively (Kresak 1982). The simultaneous comet-like appearance and 
asteroid-like orbit establish 313P as a member of the active asteroids population. 

The active asteroids are driven by a variety of mechanisms, from impact, to sublimation, 
to thermal fracture and to rotational instabilities, all previously thought to lie beyond the 
realm of observation (Jewitt 2012). In a majority of objects, the observational constraints are 
sufficiently limited that no single activity mechanism can be uniquely identihed. However, 
observations show that 313P was also active near perihelion in archival observations taken 
two orbits earlier in 2003, but not in 2004 (Hsieh et ah 2015, Hui and Jewitt 2015, Jewitt et 
al. 2015). This establishes 313P as one of only three main-belt bodies, along with 133P/Elst- 
Pizarro (Hsieh et al. 2004) and 238P/Read (Hsieh et al. 2011), to display activity recurrent 
on different orbits. Repetition in different orbits immediately rules out most of the possible 
activity mechanisms and leaves sublimation as the most likely remaining candidate (c.f. Hsieh 
and Jewitt 2006). In support of this conclusion, the distribution of surface brightness in data 
from the active epochs in both 2003 and 2014 is best matched by protracted (as opposed to 
impulsive) dust emission, as expected from comet-like sublimation. Even though a search for 
direct evidence of outgassing via the GN resonance fluorescence lines proved negative (setting 
an upper limit to the gas production rate <1.8 kg s“^; Jewitt et al. 2015), the measured 
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properties strongly suggest an origin by sublimation. 

We initiated a program of observations designed to characterize 313P and to determine 
the origin of its activity; here we present our initial results. 


2. OBSERVATIONS 

Initial observations using the 2.4 meter diameter Hubble Space Telescope (HST) were 
obtained from two orbits allocated under program GO 13864. Based on these observations, 
an additional six orbits were awarded as Director’s Discretionary Time under program GO 
14040. All measurements were taken using the WFG3 camera and the broadband F350LP 
hlter (full width at half-maximum (FWHM) = 4758A), having an effective central wavelength 
6230A when used to observe a solar-type source. The WFG3 camera consists of two charge- 
coupled devices each 2051x4096 pixels, with square pixels 0.04" on a side, giving a held of 
view 162"xl62". The image scale ranged from 42 km pixeD^ on UT 2014 October 14 to 90 
km pixeD^ on UT 2015 March 05, increasing in proportion to the geocentric distance. 

We sampled the object at roughly monthly intervals in order to follow its morpholog¬ 
ical and photometric development. Observations on UT 2015 January 07 were targeted to 
coincide with the passage of the Earth through the projected orbital plane of 313P. From 
this special viewing geometry, the perpendicular extent of the dust can be interpreted free 
from the ambiguous effects of projection. 

A preliminary discussion of the data from GO 13864 is given in Jewitt et ah (2015). 
Here, we present a full analysis of 313P from the combined HST datasets. A brief log of the 
observations is presented in Table [TJ 


3. RESULTS 
3.1. Morphology 

Figure ([^ shows a composite of the HST images at hxed angular scale (each panel is 
30" in width), with direction vectors showing the anti-solar and negative heliocentric velocity 
directions. To construct each panel, the images from each orbit of HST were shifted and 
combined to eliminate cosmic rays and other defects that are present in the individual images. 
In most cases, the cosmic ray removal was successful but residual large scale features caused 
by the parallactic smearing of held stars and galaxies sometimes remain. For example, the 
composite from UT 2015 January 07 shows a diffuse background arc that crosses the tail to 
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the east of the nucleus and which cannot be removed by digital processing. Several trailed 
galaxies are also evident in the composite from UT 2015 February 11. Such incompletely 
removed background sources limit the accuracy of surface photometry but, in most cases, 
do not detract from the conclusions to be reported here. 

The morphology of 313P changes systematically with time from showing a nucleus plus 
bright, fan-shaped tail in the early observations from 2014 October (discussed in Paper 1) 
to a thin, linear tail as the Earth passes through the orbital plane in January, towards 
an increasingly point-like appearance until the last observation on 2015 March 05 (Figure 
[^. The position angle of the tail swings counter-clockwise, following the progression of the 
projected anti-solar vector from panel to panel in the hgure. There is no evidence for a 
resolved coma about the nucleus. While the March 05 image at hrst appears stellar, very 
faint extended emission is present in the antisolar direction, showing that the nucleus was 
not completely bare. Nevertheless, it is clear from Figure ([^ that much less dust is present 
in the later images of the sequence than in the earlier ones. 

We computed syndyne/synchrone (Finson and Probstein 1968) models for each of the 
dates of observation in Table ([^. These models are shown in Figure (|^, where solid black 
lines are used to indicate syndynes (the loci of positions of particles of one size released 
from the nucleus over a range of times) and dashed black lines to indicate synchrones (loci of 
particles having a wide range of sizes released at a given time). We note that for 2015 January 
07, all synchrones and syndynes are collapsed onto the projected orbit, with synchrones later 
than 2014 September 03 to the East of the nucleus, and all earlier synchrones to the west. 
The dust particles are assumed to be released from the nucleus at zero initial speed and then 
to be accelerated by solar radiation pressure with value Pqq, where {3 is dimensionless and 
Qq is the local gravitational acceleration to the Sun. For particles whose circumference is 
larger than the wavelength of light, 27ra > A, it may be shown that (3 ~ a“^, where is 
the radius expressed in microns (Bohren and Huffman 1983). 

Figure (|^ shows that the curvature of the dust tail and the sky-plane rotation of the tail 
with respect to time are both naturally reproduced by the syndyne trajectories of particles 
having characteristic (3 ~ 0.02, or size a ~ 50 /im. The same conclusion was reached by 
Jewitt et ah (2015), Hsieh et ah (2015), Hui and Jewitt (2015) and Pozuelos et ah (2015). 
The failure of synchrone trajectories (straight lines in Figure]^ to match the tail isophotes 
is evidence that the dust emission from this body is unlikely to be caused by impact or any 
other impulsive mechanism. 
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3.2. Photometry 


We measured photometry of 313P in two ways. The brightness within a set of projected 
circular apertures of hxed angular radii is listed in Table ([^. Photometry of a distributed 
source using hxed angular radius apertures is affected by the geocentric distance, A, because 
the volume of the coma sampled by each aperture increases with the cube of increasing 
A. Therefore, we also determined the brightness within a set of apertures scaled to have 
hxed linear radius at the instantaneous distance of the object. These measurements are 
summarized in Table (|^. The sky level for all measurements was determined from the 
median count within a concentric “sky” annulus extending from 6" to 12". As discussed 
in Section (3.3) the coma surface brightness falls quickly with increasing distance from the 
nucleus, and background subtraction from the sky annulus was found to be not critical to 
the photometry. 


The brightness of 313P shows a steady decline with time. In Figure ([^ we plot the 
brightness within a set of concentric photometry apertures having hxed linear radii, from 
500 km to 6000 km, projected to the distance of 313P. The Figure shows that the apparent 
brightness of 313P faded by ~3 magnitudes between the hrst measurement on October 14 
and the last on March 05. Some of the fading is caused by the increasing distance and the 
larger phase angle of the later observations (Table [^. 

In Figure]^ we have corrected the apparent magnitudes to “absolute” magnitudes using 


Hy = V - 5 log(i?A) + 2.5 logio(*h(a)). (1) 

Here, $( 0 ;) is the phase function at phase angle a, equal to the ratio of the scattered light 
at a to that at a = 0°. We assumed the phase function formalism of Bowell et ah (1989) 
with parameter g = 0.15, as appropriate for a C-type object, as suggested by the color of 
313P and by its orbital location in the outer asteroid belt. The phase function of 313P is 
unmeasured, introducing an uncertainty into the value of Hy that is larger than the (~0.01 
magnitude) uncertainty of the photometry. To estimate the possible size of this uncertainty, 
we set the error on Hy to be equal to the difference between the phase function corrections 
for C-type and S-type objects. Over the range of phase angles at which 313P was observed, 
this difference is typically ~0.1 magnitudes. Absolute magnitudes are given in Table (|^ 
with their statistical uncertainties. 

The absolute magnitudes are related to the effective scattering cross-section of the ma¬ 
terial within the photometry aperture, Cg (km^), by 
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where pv is the geometric albedo of 313P and niQy is the apparent magnitnde of the Snn, 
both at the wavelength of the V hlter. We assnme mQy = -26.77. The resnlting scatter¬ 
ing cross-sections are listed in Table (|^, compnted assnming pv = 0.05 (c.f. Fernandez et 
al. 2013). Uncertainties on Ce are systematic in natnre, dominated by the phase fnnction 
correction as well as by the assnmption of the geometric albedo. 

The absolnte brightness of the coma annnli in Fignre Q decreases nntil the last ob¬ 
servation on March 05 (DOY 429), bnt the central apertnre flattens after UT 2015 Jannary 
21 (DOY 386) and, within the nncertainties of measnrement, remains constant. The central 
apertnre is most strongly sensitive to the brightness of the nnclens. The effective cross- 
section within the 500 km radins apertnre is = 2.2±0.1 km^ (Table [^. Formally, this 
gives an npper limit to the nnclens cross-section becanse there remains a coma even in the 
latest observations (see Fignre [^. The corresponding valne of the effective circnlar radins of 
the nnclens is = (Oe/vr)^/^ = 0.8 km. 


3.3. Radial Surface Brightness Profile 

We measnred the snrface brightness, T,{9), as a fnnction of the angnlar radius, 9, for 
each epoch of observation within a set of concentric annuli centered in the nucleus. For 
this purpose, we hrst shifted the individual images within each HST orbit to be aligned on 
the nucleus. Cosmic rays and background sources were then eliminated by computing the 
median data number within each pixel from each date. Sky background subtraction was 
determined from the median signal measured in a concentric annulus with inner and outer 
radii 9 = 6.0" and 12", respectively. Differences between the prohles, when normalized to 
unity at the peak, are modest and correspondingly difficult to illustrate graphically. We show 
the profiles on UT 2014 October 14 (Figure]^ and 2015 March 05 (Figure]^, these being 
the first and last dates of observation and illustrating the largest profile differences. The 
noise on the profiles is indicated by the scatter between measurements in adjacent annuli, 
with the profile on March 05 being very noisy as a result of the greater distance and secular 
fading of the coma. 

The surface brightness profiles all display a central bump at angular radii 9 < 0.2", 
caused by the combined effects of scattering from the nucleus and from the coma, all con¬ 
volved with the point-spread function (PSF). Our main interest in the radial surface bright¬ 
ness profiles is to estimate the effective cross-section (and size) of the nucleus. In principle. 
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the nucleus and dust contributions can be disentangled by deconvolution, using the known 
PSF. In practice, however, deconvolution acts as a noise ampliher, and we instead elected 
to estimate the nucleus contribution to the surface brightness prohle by convolving simple 
models of the intrinsic prohle with the instrumental PSF. To construct the models, we as¬ 
sumed that the intrinsic prohle can be represented by a centrally located, unresolved nucleus, 
represented by a single pixel of strength 6n{d), and a coma in which the surface brightness 
varies with angular radius 6 as S(6*) = k6~^, where A; is a constant representing the strength 
of the coma. We solved for index p by htting the observed surface brightness prohles in 
the radius range 0.2"< 6 < 1.0", this being large enough to avoid the central bump in the 
radial surface brightness. Values of p are listed together with their formal, least-squares ht 
uncertainties in Table (|^. The actual uncertainties are larger than the listed formal errors 
because, even within the tiny 0.2" to 1.0" range, deviations from power-law behavior are 
apparent. For example, the £t to data from 2015 February 11 was rendered particularly 
uncertain by structure inside 6 = 1" likely due to imperfectly removed cosmic rays and is 
omitted from further discussion. 

The model is computed from the convolution 

s(0) = psF{e) * (5„(0) + ke-p). (3) 

Here, PSF{6) is the two-dimensional representation of the PSF computed using TINYTIM 
software (Krist et ah 2011), 6n{0) is the Dirac delta function to represent the nucleus. All 
convolution calculations were done at a resolution 0.01" (pixel)”^ and then rebinned to the 
WFC3 resolution of 0.04" (pixel)”^ to compare with the data. 

Models were parametrized in terms of the quantity F, dehned by 

F = __ (4) 

K{o) + 2nkj^-^e(^-p^de 

Quantity F gives the ratio of the light scattered from the nucleus to the total light (nucleus 
plus coma) scattered within the aperture of angular radius 6 < 0.2". F = 0 corresponds to 
no nucleus contribution while F = 1 corresponds to a bare nucleus. The denominator of 
Equation (|^ is proportional to the hxed-angle photometry summarized in Table ([^. 

Results from the convolution models are summarized in Table (|^. Even models in 
which the nucleus is assumed to be negligible {F = 0) display a central surface brightness 
bump, owing to smearing of the coma brightness by convolution with the PSF. However, all 
successful hts to the data require the presence of a nucleus (i.e. F > 0) in order to match the 



height of the central bump. We obtained values of F scattered in the range 0.6 < F < 0.9, 
as listed in Table Q. The Table lists the nucleus cross-section computed from Cn = FCq, 2 , 
where Co ,2 is the cross-section within the central 0.2" taken from Table ([^. The hnal column 
of Table Q lists the effective nucleus radius in kilometers, computed from = (Cn/vr)^/^. 
The uncertainties on F are estimated from the model hts to the data. Uncertainties on Cn 
and were calculated by propagation of errors. In all cases, the uncertainties should be 
regarded as approximate. As before, all listed cross-sections are predicated on the assumption 
of geometric albedo pv = 0.05. 

Table Q shows that the effective nucleus radius varies in the range 0.6 < < 0.8 km. 

We do not know if this variability reflects uncertainties introduced by approximations in the 
method of coma removal, or by real variations associated with the rotation of an aspherical 
nucleus, or by a combination of the two. However, we note that the htted values of Cn in 
column 5 of the Table are broadly consistent with the value Cn = 2.0 km^ obtained from the 
500 km aperture photometry in UT 2015 March (Table [^. The sum of the evidence suggests 
that 0.7 km is a reasonable measure of the nucleus radius. With density p = 10^ kg m“^, 
the escape velocity from a non-rotating, spherical 0.7 km radius nucleus is W = 0.5 m s“^. 


3.4. Perpendicular Profile 

The UT 2015 January 07 composite image was used to determine the surface brightness 
prohle of the dust in the direction perpendicular to the orbit. With the Earth in the orbital 
plane on this date, the perpendicular extent accurately represents the vertical thickness of the 
dust distribution. The resulting FWHM measurements are presented in Figure ([^. There, 
vertical error bars denote ±10% photometric uncertainties in the width determination, while 
horizontal bars mark the width of the dust tail segment used to obtain the measurement, 
the latter increasing with distance from the nucleus owing to the low surface brightness of 
the dust. The dust projected to the west of the nucleus is particularly faint and therefore 
difficult to measure. As a result we show only the three FWHM measurements to the west 
in which we have conhdence. Likewise, although the dust is visible in Figure ([^ more than 
10" to the east of the nucleus, we were unable to obtain convincing FWHM measurements 
beyond this distance owing to the low dust surface brightness and interference by a smeared 
galaxy (faintly visible in the hgure as an oblique arc on the left hand side of the UT 2015 
January 07 panel in Figure [^. 

The dust trail in Figure ([^ is narrow, indicating that the dust particles are ejected at 
low velocity. For dust particles ejected by gas drag the terminal velocity is proportional to 
the inverse square-root of the particle size, V oc In this case, it can be shown that the 
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width of the resulting trail, wt-, is related to the distance from the nucleus, by 




Sir / 


( 5 ) 


where V± is the component of the ejection velocity measured in the direction perpendicular 
to the orbit and Qq is the local solar gravitational acceleration (Jewitt et ah 2014b). For 
simplicity, we assume that It is proportional to 6 , the angular distance from the nucleus 
measured in the plane of the sky. Two curves in Figure ([^ show Equation (|^ with V± = 
2.5 m s“^ (on the east side of the nucleus) and Vj_ = 5.1 m s“^ (on the west side). The 
assumption that It oc 9 neglects projection effects, and is less accurate for the much older 
particles to the west of the nucleus, as may be seen in the curved syndynes of Figure]^. As 
a result, we take V± = 2.5 m s“^ as the better estimate of the dust ejection velocity. For a 
different size-dependence of V±, the width wt becomes a function of the particle size. The 
measured trail width in this case can be taken as applicable at the optically dominant grain 
size which, from morphological analysis (section |3.1[ ), is of order 50 to 100 pm. Pozuelos et 
al. (2015) described different Finson-Probstein models in which the size dependence of the 
velocity is varied. They reported that V oc gives a better £t to the morphology than 
V oc or E oc but the effects are subtle and the signihcance of the difference 

is unclear. Regardless, unless particle ejection from the nucleus is highly collimated in the 
orbital plane, a possibility which would seem to be physically improbable, then we may 
conclude that V± is only a few meters per second, and that it provides a measure of the total 
velocity of ejection of the dust. 

The low dust velocity accounts for the non-detection of coma in 313P. A particle ejected 
towards the Sun at speed V has a turn-around distance 


2 / 3 ^©' 


( 6 ) 


Substituting V = 2.5 m s“^, fi = 0.01 to 0.02 and Qq = 10“^ m s“^ gives s = 160 to 310 
km, corresponding to only 0.1" to 0.2" at 2 AU. While technically resolvable with HST, 
this near-nucleus region of the surface brightness prohle is dominated by the point-spread 
function of the vastly brighter nucleus (e.g. see Figures]^ and [^, making any coma difficult 
to detect. 


We conclude that the dust is ejected from 313P at a speed two orders of magnitude 
smaller than the sound speed (~400 m s“^) in sublimated gas at ~3 AU, but slightly larger 
than the likely R. = 0.5 m s“^ gravitational escape speed from the nucleus. Very small 
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launch speeds have been found in other active asteroids, most notably in the probable ice- 
sublimator 133P/Elst-Pizarro (Hsieh et ah 2004, Jewitt et ah 2014b). Low terminal velocities 
can result from weak gas flow, as expected for ice at distances ~ 3 AU. Low terminal 
velocities can also be produced by the small size of the sublimating ice patch, which reduces 
the acceleration length over which gas drag can act (Jewitt et al. 2014b). 


4. DISCUSSION 

From Table ([^, it may be seen that the amount of dust in 313P (represented by Cd = 
(76(6000) — Cn) decreases from October to March by a factor of ~5. During this time, the 
heliocentric distance increased by only 8%, from 2.4 to 2.6 AU (Table and the sub-solar 
equilibrium temperature (which scales as fell by only ~4%. The implied distance 

dependence (dust cross-section oc is unreasonably steep, and we conclude that the 

fading is not caused by the changing distance. The data in Figure (|^ show that the coma 
absolute magnitudes have e-folding timescales ~75 days, regardless of the radius of the 
photometry annulus used to measure the coma. What determines the long timescale of the 
fading of 313P? 


4.1. Radiation Pressure Sweeping 

Radiation pressure is capable of sweeping particles from the vicinity of the nucleus, 
but a simple calculation shows that the sweeping timescales are short compared to the 
~75 day fading timescale. To see this, we note that the dominant particles as judged from 
syndyne/synchrone analyses have radiation pressure acceleration parameter (3 ~ 0.01 to 0.02 
(corresponding to particles of ~50 to 100 /xm size). Neglecting the initial velocity of ejection, 
the time taken by solar radiation pressure to accelerate a dust particle over a distance L, is 
given by 


r = 



(7) 


where Qq is the local solar gravitational acceleration. At 2.5 AU from the Sun, the solar 
gravitational acceleration is Qq = 10“^ m s“^. We consider, for example, dust in the 6000 
km radius photometry aperture. With L = 6x10® m. Equation ([^ gives r ~ 8 to 11 xlO® s 
(9 to 13 days) and these are upper limits because of the neglect of the initial velocity. This 
rapid clearing shows that radiation pressure is unlikely to produce the decline in the coma 
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cross-sections that occurs on timescales of ~75 days. Indeed, the persistence of the coma 
over hve months requires that dust be continually released from the nucleus, evidently at 
a declining rate. This photometry-based conclusion is consistent with the inference, made 
using Finson-Probstein dust dynamics models, that the mass loss occurred over a protracted 
period (Jewitt et ah 2015, Hsieh et ah 2015, Hui and Jewitt 2015, Pozuelos et ah 2015) and 
with an origin by sublimation. 

Given that the dust is replenished on the timescale given by Equation ([^, we estimate 
the dust mass loss rate from the coma using 


dMd _ paCd 

dt T 


( 8 ) 


where p = 10^ kg m“^ is again the assumed bulk density of ejected grains, a is their effective 
radius, Cd is the scattering cross section in the coma (Table and r is from Equation 
([^. We use a = 50 to 100 /im, as found from the syndyne/synchrone models. Substitution 
into Equation ([^ gives values dMd/dt = 0.6 to 0.9 kg s“^ in 2014 October declining to 
~0.2 to 0.3 kg s“^ by 2015 March. These estimates compare with an upper limit to the gas 
production rate < 1.8 kg s“^ based on the non-detection of emission lines in a spectrum taken 
on UT 2014 October 22 (Jewitt et al. 2015). Peak dust production rates were independently 
estimated by Pozuelos et al. (2015) at dMd/dt = 0.2 to 0.8 kg s“^ (UT 2014 September 21). 
Given the many uncertainties in the dust and model parameters, we consider these values in 
good agreement. 


4.2. Pit Source 

Instead, the long timescale of the coma fading leads us to consider illumination effects, of 
which we distinguish three types. First, the orbit of 313P is eccentric, causing the heliocentric 
distance to vary from 2.391 AU at perihelion to 3.921 AU at aphelion. This gives a factor 
~2.7 in solar insolation, which might be expected to have a measurable impact on the rates 
of sublimation at perihelion vs. aphelion. However, Table ([^ shows that the change in the 
heliocentric distance in any 75 day interval is negligible and, therefore, a decrease in the 
insolation is unlikely to be responsible for the observed fading. Second, a seasonal variation 
of the solar insolation at any point on the surface of 313P will result from non-zero obliquity 
of the spin, being largest for obliquity = 90°, when the Sun can cross both the projected 
equator and the poles of 313P at different points in the orbit. Third, the pattern of local 
shadows on the surface will be modulated by the position of 313P in its orbit. Since the 
sublimating area on 313P is very small (Jewitt et al. 2015) there is reason to expect that 
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local shadows may play an important role in modulating the mass loss. 

We consider a simple model to attempt to capture the second and third of these effects. 
In this model, the instantaneous equilibrium sublimation rate, Fg (kg m“^ is calculated 
from the energy balance equation 


cos(^(f)) = eaT^ + L(T)F,(f) (9) 

rjj 

in which Fo = 1360 W m ^ is the Solar constant, A is the Bond albedo, e is the emissivity, 
a is the Stefan-Boltzmann constant, L{T) is the latent heat of sublimation for water ice at 
temperature, T, and Fg is the equilibrium sublimation flux (kg m“^ s“^). We took A = 0.05, 
e = 0.9, a = 5.67x10“* W m“^ K“'^, while L{T) was obtained from Washburn (1926) and 
all calculations were performed for th = 2.5 AU, representative of the heliocentric distance 
of 313P (Table [^. 

The angle i{t) in Equation ([^ is the zenith angle of the Sun as viewed from the subli¬ 
mating surface. Angle i (which must be < 90° for Equation ([^ to hold) is determined by the 
instantaneous orientation of the nucleus spin vector relative to the Sun direction, by ifj, the 
latitude of the sublimating patch, by the local slope of the surface and by the time of day. 
Neither the shape nor the spin vector of 313P are well known. In the following, we assume 
that the obliquity is 90° in order to maximize the seasonal effects. In this geometry, and 
assuming a spherical nucleus with a spin vector perpendicular to the Sun-nucleus line, the 
local solar zenith angle may be written as cos{i(t)) = cos('0) cos(ci;t), where u = 27i/P is the 
angular frequency of the rotation, P is the rotation period and t is the time. We arbitrarily 
take P = 4 hr, although the period is unimportant for the relative effects described here. 

Solutions to Equation (|^ are plotted in Figure (|^ for the diurnal variation of the ice 
sublimation rate on the rotational equator {'ip = 0°) and at '0 = 45°, as solid red and blue 
curves, respectively. The hgure shows that the largest sublimation rates at 2.5 AU are 
Fg ^ 5 X 10“® kg m“^ s“^ on the equator, falling to about 60% of this value at 45° latitude. 
As expected, sublimation at 45° is limited in duration as well as in rate relative to the 
equatorial case. Continued mass loss from a region of the surface at the rate Fg (kg m“^ 
s“^) will cause the sublimating surface to sink relative to the surrounding non-sublimating 
surface at rate di/dt = Fg/p, where p is the bulk density of the surface layers. With p = 10* 
kg m“*, the surface recession rate is d£/dt = 5 x 10“* m s“*, or 4 mm day“*. Progressive 
sublimation will create a pit whose floor remains partly self-shadowed, leading to variations 
of the sublimation rate on both diurnal and seasonal timescales. Sublimation pits have 
long been discussed in the context of active comets (e.g. Keller et al. 1994), although their 
formation mechanism is uncertain (Thomas et al. 2013). They have recently been imaged in 
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spectacular detail on the nucleus of comet 67P/Churyumov-Gerasimenko (Figure [^). 

We explored a simple model of self-shadowing, in which the pit is rectangular in cross- 
section with a depth-to-diameter ratio d/D (see Figure]^). We assume that sublimation 
proceeds from the floor of the pit, although there is some suggestion that at least one 
pit on 67P sublimates through its walls (Vincent et al. 2015). The detailed geometry is 
not crucial to the argument that follows; what matters most is that heating of the ice is 
topographically obstructed, and this depends mainly on the ratio of the depth (or height) 
to the horizontal scale of the sublimating patch. For the assumed pit source, the fraction of 
the floor illuminated by the Sun is 


/(t) ~ 1 - 


d 

D tan(i(t)) 


( 10 ) 


provided i > ic 
compute 


tan ^{d/D), and f{t) = 0 otherwise. Equations ([^ and (10) were used to 


AM = r f{t)F,{t)dt (11) 

Jo 


where the integral is taken over one rotation of the nucleus. Equation ( pT| gives the total 
mass lost in one nucleus rotation per unit area by equilibrium sublimation from a surface 
element at latitude 'ip. 


Figure ( p!o| shows a comparison between the photometry and the model results from 
Equation ( pT| . For clarity of presentation we show only measurements from the 500:1000 
km annulus, but other annular photometry measurements give consistent results, since they 
show similar rates of fading. To make the hgure we have assumed that the sublimated mass 
from Equation ( pT| and the scattering cross-section (and hence the measured brightness) 
are proportional, AM oc C^- We hnd that models with d/D = Q (i.e. surface sublimation) 
and d/D = 1 are strongly inconsistent with the measured fading rates because they predict 
brightness variations with mean anomaly that are, respectively, too flat and too steep to 
match the observations. On the other hand, models with d/D ~ 1/2 better match the rate 
of decline in the brightness of 313P by virtue of the inclusion of self-shadowing in the pit. 


The curves in Figure (10) are non-unique, and the model behind them is highly sim- 
plihed. For instance, we are forced to assume the obliquity in the absence of constraining 
evidence and we assumed that the nucleus of 313P is spherical when such a small body is un¬ 
likely to be so. A different obliquity and a different nucleus body-shape would give different 
results. We have also neglected any possibility of an insulating dust mantle on the sublimat- 
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ing ice at the bottom of the pit. This is likely a serious omission, in that pit sublimation may 
be stifled both by self-shadowing and by the accumulation of debris (see Guilbert-Lepoutre 
et al. 2015 for a modern discussion of these effects). Nevertheless, our model is sufficient to 
show that seasonal effects caused by non-zero obliquity and local shadowing in a pit source 
can easily match the observed fading rate, and give a physically plausible pit geometry. 
This conclusion can be extended to other active asteroids, for instance 133P/Elst-Pizarro 
(Hsieh et al. 2004), in which activity is known to be restricted to a substantial range of true 
anomalies. 

For an equatorially located pit with d/D = 1/2, Figure (|^ indicates a rotationally- 
averaged mass loss rate Fg ~ 1x10“® kg m“^ s“^. To compare with the estimates of the 
dust mass production rate (near 1 kg s“^ at peak), we need to know the ratio of dust to 
gas production rates, fdg- Measurements of active comets generally give values fdg > 1 (this 
is possible because the escaping gas travels much faster than the ejected dust, maintaining 
momentum equipartition between the two components). For example, observations of comet 
2P/Encke give 10 < fdg < 30 (Reach et al. 2000). The area of exposed ice needed to supply 
a dust mass loss rate dMd/dt is 


nrl 


1 \ dMd 

fdgFs) dt 


( 12 ) 


With dMd/dt = 1 kg s“^, as inferred above, Fg = 1x10“® kg m“^ s“^ and conservatively 
taking fdg = 10, we find vrr^ = 10^ m^, and Vg = 56 m. This sublimating area corresponds 
to ~0.2% of the surface of a 700 m radius spherical nucleus. It is possible that the source is 
a single region of radius and depth d = ~ 56 m, perhaps formed by a small impact. Such 

a pit would take a time dp/{Fg) ~ 6x10® s (200 years) to grow. We think it more likely 
that the sublimating area consists of a number of pits, each less deep and younger than this 
estimate. For example, N equal-size pits of radius and depth r^ 56 /N^F would satisfy the 
observations. With N = 100, each pit would be ~6 m deep and have an excavation time 
~20 years. Unfortunately, we possess no observational constraint on N. 

Low dust velocities in 133P were interpreted as evidence for sublimation from patches of 
limited size, because the acceleration length and the terminal velocity for gas-entrained dust 
particles scale with the physical size of the sublimating source region (Jewitt et al. 2014b). 
Dust in 313P is ejected with small velocities comparable to those in 133P (c.f. Figure]^, and 
sublimation from a set of small pits may again be responsible. It may be natural to expect 
that an exposed ice surface should develop into a set of pits in response to fallback mantling, 
the presence of surface irregularities such as boulders and recession of the sublimating surface 
beneath the physical surface of the adjacent nucleus. 




15 


4.3. Variation Around the Orbit 

Seasonal effects in all three of the active asteroids that have displayed repetitive activity 
in different orbits are plotted in Fignre ( [IT| ). These are the strongest candidates for being 
driven by the snblimation of ice. The data were compiled for 133P from Hsieh et ah (2010), 
for 238P from Hsieh et ah (2011) and for 313P from Jewitt et ah (2015), Hsieh et ah (2015), 
Hni and Jewitt (2015) and from the present work. There is an absence of activity within 
the trne anomaly range v ~ 180±90° for each object. Activity is observed to abont ±60° in 
238P. Prototype 133P shows activity mainly after perihelion, np to i/ ~ 80°, snggesting the 
action of a thermal lag (Hsieh et ah 2010). The variation of activity aronnd the orbit shown 
in Fignre ( [IT| ) is consistent with seasonal modnlation on all three objects, as expected if 
snblimating ice is responsible. However, it shonld be pointed ont that the fnll range of trne 
anomaly is incompletely sampled, especially for 238P and 313P, and it is possible that activity 
occnrs over a wider fraction of each orbit than existing observations reveal. In addition, it is 
likely that observational bias favors the detection of objects near perihelion becanse objects 
are brighter there than at aphelion. For example, typical perihelion and aphelion distances 
of the three above objects are ~2.5 AU and 4.0 AU, respectively. Observed as point sonrces 
at opposition (corresponding to geocentric distances ~1.5 AU and 3.0 AU, respectively), 
the inverse sqnare law predicts the objects to be fainter at aphelion than perihelion by 
~2.5 magnitndes. In addition, the ratio of the eqnilibrinm snblimation rates at perihelion 
and aphelion is an order of magnitnde or more, constitnting an observationally formidable 
obstacle to the detection of activity at large distances. More work is needed to buttress the 
apparent concentration of activity in these bodies near perihelion. 


4.4. Timescales 

The inferred mass loss rate dMd/dt ~ 1 kg s“^, if sustained over 5 months, would 
correspond to a total mass lost from 313P of AMd ~10^ kg per Pk = 5.6 year orbit. For 
comparison, the mass of the nucleus, represented as a sphere of radius = 700 m and 
density p = 10^ kg m“^ is ~ 10^^ kg. Thus, the fractional mass loss per orbit is 
AMd/M^ 10 ^ and 313P could sustain continuous activity at the present rate for a time 
PKMn/AMd ~ 6x10^ yr. The age of 313P is unknown but is likely to be much greater. 
Hsieh et al. (2015) report a dynamical association with the 160±35 Myr old Lixiahua asteroid 
family (c.f. Novakovic et al. 2010) which, if real, would suggest that 313P is about 100 times 
older than its sublimation age. The approximate collisional destruction time for a = 700 
m radius asteroid is even older at 10® yr (Bottke et al. 2005). The survival of ice in 313P for 
10® and 10® yr periods can be simply explained by the presence of a refractory mantle thick 
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enough to stifle sublimation of near-surface ice (Schorghofer 2008). Such a mantle would 
block sublimation until penetrated, perhaps by impact or by surface instability (Hsieh and 
Jewitt 2006). The duty cycle for activity (i.e. the fraction of the time for which the body is 
active) need only be < 10“^ for ice to survive for the age of the Lixiahua family or < 10“^ 
to survive for the collisional lifetime of the object, respectively. A duty cycle of < 10“^ or 
< 10“^ implies that, for every observed active case like that of 313P, there are > 10^ or 
> 10^ similar but dormant, ice-containing bodies in the asteroid belt. 
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5. SUMMARY 

313P/Gibbs is an active asteroid ejecting dnst near perihelion in two different orbits, 
suggesting an origin by the intermittent sublimation of near-surface ice. We have used the 
Hubble Space Telescope to study this object in detail, in order to better constrain the nature 
of its activity. We find that 

1. The nucleus radius, estimated from convolution models of the surface brightness prohle, 
lies in the range 0.6 to 0.8 km (geometric albedo 0.05 assumed). 

2. The dust distribution is consistent with the continued ejection of large particles (radi¬ 
ation pressure parameter (3 = 0.01 to 0.02, corresponding to particle radii ~50 to 100 
fim) in observations extending from 2014 October to 2015 March. 

3. In-plane observations show that these particles are ejected slowly, with characteristic 
velocities normal to the orbit plane of ~2.5 m s“^. 

4. The peak mass loss rate in dust is of order 1 kg s“^, decreasing on an e-folding time 
~75 day. The area of exposed ice needed to supply this rate is ~10"^ m^, corresponding 
to only ~0.2% of the nucleus surface. 

5. The secular fading is too slow to be caused by radiation pressure sweeping of dust 
particles from the coma but too fast to be related to the (marginally) increasing he¬ 
liocentric distance. Instead, we show by a simple model that the fading timescale is 
consistent with protracted sublimation of ice from a topographically shadowed region, 
for example in a pit having depth-to-diameter ratio ~0.5. Such pits are a natural 
product of localized ice sublimation. 

6. The ejected dust mass is ~10^ kg per orbit, approximately 10“^ of the 10^^ kg nucleus 
mass. If comparably active in every orbit, the mass-loss lifetime is conhned to ~0.6 
Myr, far smaller than the ~160±35 Myr age of the Lixiahua family of which 313P is a 
likely member, and smaller than the ~1 Gyr lifetime to collisions. To reconcile these 
timescales requires that the duty cycle for sublimation-driven mass-loss be ^10“^ to 
<10“^, respectively. 

Based in part on observations made with the NASA/ESA Hubble Space Telescope, with 
data obtained via the Space Telescope Science Institute (STSGI). Support for program 14040 
was provided by NASA through a grant from STSGI, operated by AURA, Inc., under contract 
NAS 5-26555. We thank Linda Dressel, Alison Vick and other members of the STScI ground 
system team for their expert help and the anonymous referee for comments. 
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Table 1. Observing Geometry 


UT Date and Time 

DOY*^ 

ATpb 

z/" 


A" 





2014 Oct 14 13:12 - 13:48 

287 

47 

14 

2.405 

1.451 

8.9 

10.5 

247.9 

7.5 

2014 Oct 28 21:25 - 23:16 

301 

61 

18 

2.415 

1.522 

13.0 

37.1 

248.8 

6.8 

2014 Dec 01 04:36 - 05:13 

335 

95 

28 

2.446 

1.835 

21.0 

58.8 

249.2 

3.6 

2014 Dec 15 07:36 - 08:12 

349 

109 

32 

2.463 

2.011 

22.6 

62.9 

248.7 

2.1 

2015 Jan 07 16:20 - 18:08 

372 

132 

38 

2.495 

2.326 

23.2 

67.7 

247.6 

-0.0 

2015 Jan 21 03:07 - 03:43 

386 

146 

42 

2.516 

2.513 

22.6 

69.9 

247.1 

-1.0 

2015 Feb 11 19:37 - 20:14 

407 

167 

47 

2.553 

2.809 

20.5 

73.3 

246.8 

-2.2 

2015 Mar 05 14:00 - 15:48 

429 

189 

53 

2.593 

3.085 

17.5 

76.9 

247.1 

-3.0 


'^Day of Year, UT 2014 January 01 = 1 

'^Number of days past perihelion (UT 2014 August 28 = DOY 240) 

'^True anomaly, in degrees 
'^Heliocentric distance, in AU 
^Geocentric distance, in AU 
^Phase angle, in degrees 

^Position angle of the projected anti-Solar direction, in degrees 

'^Position angle of the projected negative heliocentric velocity vector, in degrees 

'Angle of Earth above the orbital plane, in degrees 
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Table 2. Photometry with Fixed Angle Apertures 


UT Date 


niY^ 

Hv^ 

Ce[km‘^]^ 

Oct 14 

0.2 

20.62±0.01 

17.3±0.07 

3.3±0.2 

Oct 14 

1.0 

20.02±0.01 

16.7±0.07 

5.8±0.4 

Oct 14 

4.0 

19.60±0.01 

16.3±0.07 

8.5±0.6 

Oct 14 

6.0 

19.40±0.01 

16.1±0.07 

10.2±0.7 

Oct 28 

0.2 

20.97±0.01 

17.4±0.09 

3.1±0.3 

Oct 28 

1.0 

20.36±0.01 

16.8±0.09 

5.4±0.5 

Oct 28 

4.0 

19.81±0.01 

16.2±0.09 

9.0±0.8 

Oct 28 

6.0 

19.72±0.01 

16.1±0.09 

9.8±0.9 

Dec 01 

0.2 

21.83±0.01 

17.54±0.12 

2.7±0.3 

Dec 01 

1.0 

21.27±0.01 

16.98±0.12 

4.5±0.5 

Dec 01 

4.0 

20.63±0.01 

16.34±0.12 

8.0±1.0 

Dec 01 

6.0 

20.43±0.01 

16.14±0.12 

9.7±1.2 

Dec 15 

0.2 

22.17±0.01 

17.62±0.13 

2.5±0.3 

Dec 15 

1.0 

21.66±0.01 

17.11±0.13 

3.9±0.5 

Dec 15 

4.0 

21.01±0.01 

16.46±0.13 

7.2±0.9 

Dec 15 

6.0 

20.85±0.01 

16.30±0.13 

8.3±1.1 

Jan 07 

0.2 

22.67±0.01 

18.12±0.13 

1.6±0.2 

Jan 07 

1.0 

22.08±0.01 

17.53±0.13 

2.7±0.3 

Jan 07 

4.0 

21.28±0.01 

16.73±0.13 

5.6±0.7 

Jan 07 

6.0 

21.01±0.01 

16.46±0.13 

7.2±0.9 

Jan 21 

0.2 

22.96±0.01 

17.87±0.13 

2.0±0.3 

Jan 21 

1.0 

22.48±0.01 

17.39±0.13 

3.1±0.4 

Jan 21 

4.0 

21.78±0.01 

16.69±0.13 

5.8±0.8 

Jan 21 

6.0 

21.74±0.01 

16.65±0.13 

6.0±0.8 

Feb 11 

0.2 

23.05±0.01 

17.76±0.13 

2.2±0.3 

Feb 11 

1.0 

22.68±0.01 

17.39±0.13 

3.1±0.4 

Feb 11 

4.0 

22.18±0.01 

16.89±0.13 

4.8±0.6 

Feb 11 

6.0 

22.06±0.01 

16.77±0.13 

5.4±0.7 

Mar 05 

0.2 

23.27±0.01 

17.83±0.11 

2.0±0.2 

Mar 05 

1.0 

22.57±0.01 

17.53±0.11 

2.7±0.2 

Mar 05 

4.0 

22.63±0.01 

17.19±0.11 

3.7±0.4 

March 05 

6.0 

22.60±0.01 

17.16±0.11 

3.8±0.4 


^Angular radius of photometry aperture, in arcsec 
^Apparent V magnitude 

‘^Absolute magnitude computed from Equation 

^Effective scattering cross-section computed from Equation 
1 ^, km^ 






Table 3. 


Photometry with Fixed Linear Apertures 


UT Date 

Quantity^ 

500 km 

1000 km 

2000 km 

4000 km 

6000 km 

Oct 14 

V 

20.30 

20.04 

19.77 

19.53 

19.44 

Oct 14 

Hv 

16.99 

16.73 

16.46 

16.22 

16.13 

Oct 14 

Ce 

4.4 

5.6 

7.2 

9.0 

9.7 

Oct 28 

V 

20.67 

20.39 

20.11 

19.77 

19.67 

Oct 28 

Hv 

17.08 

16.80 

16.52 

16.18 

16.08 

Oct 28 

Ce 

4.1 

5.3 

6.8 

9.3 

10.2 

Dec 01 

V 

21.64 

21.38 

21.10 

20.79 

20.59 

Dec 01 

Hv 

17.35 

17.09 

16.81 

16.50 

16.30 

Dec 01 

Ce 

3.2 

4.0 

5.2 

6.9 

8.3 

Dec 15 

V 

21.99 

21.79 

21.55 

21.25 

21.05 

Dec 15 

Hv 

17.44 

17.24 

17.00 

16.70 

16.50 

Dec 15 

Ce 

2.9 

3.5 

4.4 

5.8 

6.9 

Jan 07 

V 

22.57 

22.31 

22.03 

21.67 

21.45 

Jan 07 

Hv 

17.65 

17.39 

17.11 

16.75 

16.53 

Jan 07 

Ce 

2.4 

3.1 

4.0 

5.5 

6.7 

Jan 21 

V 

22.86 

22.67 

22.46 

22.10 

21.94 

Jan 21 

Hv 

17.78 

17.59 

17.38 

17.11 

16.86 

Jan 21 

Ce 

2.1 

2.5 

3.1 

4.0 

5.0 

Feb 11 

V 

23.01 

22.85 

22.70 

22.52 

22.40 

Feb 11 

Hv 

17.71 

17.55 

17.40 

17.22 

17.10 

Feb 11 

Ce 

2.3 

2.6 

3.0 

3.6 

4.0 

Mar 05 

V 

23.24 

23.13 

23.01 

22.85 

22.79 

Mar 05 

Hv 

17.80 

17.69 

17.57 

17.41 

17.35 

Mar 05 

Ce 

2.1 

2.3 

2.6 

3.0 

3.2 


= apparent V magnitude, Hy = absolute V magnitude, Ce = effective scattering 
cross-section in km^ 



Table 4. Surface Brightness Fits 


UT Date 

pa 

py b 


d 

rp ® 

' n 

Oct 14 

1.60T0.01 

3.3±0.3 

0.66±0.10 

2.2T0.4 

0.8±0.1 

Oct 28 

1.60T0.01 

3.1±0.3 

0.66±0.10 

2.0T0.4 

0.8±0.1 

Dec 01 

1.52T0.01 

2.7±0.3 

0.80±0.10 

2.2T0.4 

0.8T0.1 

Dec 15 

1.77±0.01 

2.5±0.3 

0.66±0.10 

1.7±0.3 

0.7T0.1 

Jan 07 

1.42T0.01 

1.6±0.2 

0.70±0.10 

1.1T0.2 

0.6T0.1 

Jan 21 

1.65±0.01 

2.0±0.3 

0.64±0.10 

1.3T0.3 

0.6T0.1 

Feb ll’ 

~ 

- 

- 

~ 

~ 

Mar 05 

2.08T0.10 

2.2±0.2 

0.80±0.10 

1.8T0.3 

0.7±0.1 


'^Surface brightness gradient index measured in the angular radius 
range 0.2" < 6 < 1.0" 

’^Cross-section (km^) inside a 0.2" radius circle, from Table 

'^Fraction of Co ,2 contributed by the nucleus, determined from the 
convolution model 

'’Nucleus cross-section (km^) 

^Effective nucleus radius, = (Cn/7r)’/^, (km) 

’No reliable £t for p was possible using data from this date 
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Fig. 1.— Composite of HST images of 313P/Gibbs in which each panel has North to the 
top, East to the left, and is 30" wide. The dates of the images are indicated. Arrows show 
the directions of the projected antisolar vector (-S) and the negative projected heliocentric 
velocity vector (-V). Nnmbers in the lower right corner of each panel show the angle of the 
Earth above the orbital plane of the object, in degrees. 
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Right Ascension Offset [arcsec] 


Fig. 2.— Syndynes (solid black lines) for /3 = 0.1, 0.05, 0.02, 0.01, 0.005, 0.002, 0.001, 0.0005, 
0.0002 and 0.0001 are shown. However, only the (3 = 0.01 and 0.001 syndynes are labeled 
for clarity. Labels for synchrones (dashed straight lines) are also suppressed for clarity but 
can be determined as follows. The blue diamond corresponds to ejection on 2014-July-9, 
with sychrones plotted every 10 days afterwards (anti-clockwise) and 50 days before. The 
red diamond is for ejection on 2014-Aug-28 and plotting intervals are the same. 
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Day of Year 2014, DOY 


Fig. 3.— Apparent magnitudes determined within circular apertures of fixed outer radii, 
from 500 km to 6000 km, as marked. The statistical measurement error bars (±0.01 magni¬ 
tude) are too small to be seen at the scale of the plot. 
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Day of Year 2014, DOY 


Fig. 4.— Absolute magnitudes determined within circular annuli scaled to have hxed linear 
inner and outer radii, as indicated, in kilometers. The error bars on the smallest annulus, 
77^/(0:500), show systematic uncertainties resulting from the unmeasured phase function 
of 313P. Similar uncertainties exist for the larger annuli but are not shown to avoid plot 
confusion. 










Fig. 5.— Normalized surface brightness profiles of 313P measured on UT 2014 October 14 
(green circles) compared with models. Quantity F is the fraction of the signal within angular 
radius 6 < 0.2" that is contributed by scattering from the nucleus (see Equation]^. 











Normalized Surface Brightness 
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Fig. 6. 


Same as Figure (|^ but for observations taken UT 2015 March 05. 
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Fig. 7.— Full width at half maximum of the dust emitted from 313P observed on UT 2015 
January 07 from a vantage point in the orbital plane of the object. Lines to the east mark 
ejection velocities 2.5±0.5 m s“^ and to the west 5.1±0.5 m s“^. 
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Fig. 8.— a) Oblique image of a deep sublimation pit (depth to diameter ratio ~1) on 
the nucleus of Jupiter-family comet 67P/Churyumov-Gerasimenko recorded from the ESA 
Rosetta spacecraft (credit: ESA/Rosetta/MPS/OSIRIS) b) sketch dehning parameters of 
the pit model in Equation (11). 
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Sun Azimuth Angle [degree] 


Fig. 9.— Sublimation from ice exposed at two latitudes (-0 = 0° and 45°) and two pit depths 
{d/D = 0 and 0.5), as marked, as a function of the instantaneous solar azimuth angle. 



33 



True Anomaly [degree] 


Fig. 10.— The absolute magnitude, ify(500:1000) (green circles) as a function of the true 
anomaly. Three pit source sublimation models having depth to diameter ratios d/D = 
0.0 (red dash-dot line), 0.5 (black solid line) and 1.0 (blue dashed line) show the effect of 
self-shadowing for an equatorial source on a nucleus with 90° obliquity. 
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Fig. 11.— Activity state of the three repetitively active asteroids as a function of the true 
anomaly at the time of observation. Detections of activity are marked using hlled squares 
while non-detections are marked with unhlled squares. 



















